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FUNDAMENTAL PROPERTIES OF STARS USING ASTEROSEISMOLOGY FROM KEPLER & COROT AND 

INTERFEROMETRY FROM THE CHARA ARRAY 
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ABSTRACT 

We present results of a long-baseline interfcromctry campaign using the PAVO beam combiner at the 
CHARA Array to measure the angular sizes of five main-sequence stars, one subgiant and four red giant 
stars for which solar-like oscillations have been detected by either Kepler or CoRoT. By combining 
interferometric angular diameters, Hipparcos parallaxes, asteroseismic densities, bolometric fluxes and 
high-resolution spectroscopy we derive a full set of near model-independent fundamental properties 
for the sample. We first use these properties to test asteroseismic scaling relations for the frequency 
of maximum power (f ma x) and the large frequency separation (A^). We find excellent agreement 
within the observational uncertainties, and empirically show that simple estimates of asteroseismic 
radii for main-sequence stars are accurate to < 4%. We furthermore find good agreement of our 
measured effective temperatures with spectroscopic and photometric estimates with mean deviations 
for stars between T off = 4600 - 6200 K of -22 ± 32 K (with a scatter of 97 K) and -58 ± 31 K (with 
a scatter of 93 K), respectively. Finally we present a first comparison with evolutionary models, 
and find differences between observed and theoretical properties for the metal-rich main-sequence 
star HD 173701. We conclude that the constraints presented in this study will have strong potential 
for testing stellar model physics, in particular when combined with detailed modelling of individual 
oscillation frequencies. 

Subject headings: stars: oscillations — stars: late-type — techniques: photometric — techniques: 
interferometric 
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1. INTRODUCTION 

The knowledge of fundamental properties such as tem- 
perature, radius and mass of stars in different evolu- 
tionary phases plays a key role in many applications 
of modern astrophysics. Examples include the im- 
provement of model physics of stella r structure and 
evolu t ion such as convectio n (see, e.g. iDemarque et all 
19861: IMonteiro et all llIM IDeheuvels fc Michell l201lT 



Trampedach fc Steinl l201lT iPiau et al l 120111 ). the cal 
ibration of empirical relations such a s the color- 
temperature scale for cool stars (see, e.g.. lFlowerl[l996l : 
iRamfrez fc Melendezll2005l: iCasagrande et al.ll2010l ). and 
the characterization of physical pr operties and habit- 
able zones of exoplanets (see, e.g., [Baines et al.l 120081 : 
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Ivan Belle fc von Braunll200l Ivon Braun et al.ll2011al lH). 

Many methods to determine properties of single field 
stars are indirect, and therefore of limited use for im- 
proving stellar models. Asteroseismology, the study 
of stellar oscillations, is a powerful method to deter- 
mine properties of solar-type stars such as the mean 
stellar density with little model dependence (see, e.g. , 
Brown fc Giliilandl H99l lChristensen-"P alsgaard l200l 



Aerts et al.l l2010f) 7 Additionally, long-baseline interfer- 
ometry can be used to measure the angular sizes of 
stars which, in combination with a parallax, yields a 
linear radius and, when combined with an estimate for 
the bolometric flux, provides a direct measurement of a 
star's effective temp e rature (see, e.g.. iCode et al.l [19761: 
Bovaiian et al.|[2009t iBaines et al.ll2009t iBovaiian et al.l 



2012 



man et al. zuuy; mine s et al.ll/UUbit ITJovauan et al 
al|bl: iCreevev et al.l I2012T ). Therefore, the combina- 



tion of both methods in principle allows a determina- 
tion of radii, masses and temperatures of stars with little 
model dependence. 

While the potential of combining asteroseismology 
and interferometry has been long recognized (see, e.g, 
ICunha et al.l [2007), observational constraints have so 
far restricted an appli cation for cool stars to rela- 
tively few bright objects ([North et al.ll2007t iBruntt et al.l 
l2010tlBazot et al.1l2011l ). Recent technological advances, 
however, have changed this picture. The launches of 
the space telescopes CoRoT ( Convection , Rotation and 
planetary Transits. IBaglin et al.l [2006atf a) and Kepler 
([Borucki et all I201CH : iKoch et al.l l20ldh has increased 
the number of stars with detected solar-like oscillations 
to several thousands, providing a large sample span- 
ning from t he main-sequence to He-core burning red 
giant stars (iPe Ridder et al.l l2009t jHekker et al.l 120091: 
iGilliland et al.ll2010bl : iChaplin et atlboill l At the same 
time, the development of highly s ensitive instruments 
such as the PAVO be am combiner (llreland et al.l I2008T ) 
at the CHARA Array (jten Brummelaar et al.ll2005l ) have 
pushed the sensitivity limits of long-baseline interferom- 
etry, bringing into reach the brightest objects for which 
high-quality space-based asteroseismic data are available. 
Using these recent advances, we present a systematic 
combined asteroseismic and interferometric study of low- 
mass stars spanning from the main-sequence to the red 
clump. 

2. TARGET SAMPLE 

Our target sample was selected to optimize the combi- 
nation of asteroseismology and interferometry given the 
observational constraints, while also covering a large pa- 
rameter space in stellar evolution. The majority of our 
stars were taken from the sample analyzed by the Kepler 
Asteroseismic Science Consortium (KASC). We selected 
four unevolved stars, which are among the brightest oscil- 
lating solar-type stars observed b y Kepler. Note tha t our 
interferometr ic results for flCyg (Guzik et al.ll201lT ) and 
16 Cyg A&B ([Metcalfe et al.ll2012l ) will be presented else- 
where. For the Kepler giant sample, four of the brightest 
red giants with the best Hipparcos parallaxes were se- 
lected. Finally, the main-sequence stars HD 175726 and 
HD 181420 in our sample are located in the CoRoT field 
towards the galactic center, and were among the first 
CoRoT main-sequence stars with dete cted oscillations 
(|Barban et al.l 120091: iMosser et al.ll2009D . Note that our 
PAVO campaign is also targeting solar-type oscillators in 
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Fig. 1. — H-R diagram with the position of the target stars cal- 
culated using spectroscopy, photometry and Hipparcos parallaxes. 
Solar metallicity BaSTI evolutionary tracks from 0.6-2.6 Mq in 
steps of 0.1 Mq are shown as grey lines. The dashed line marks 
the approximate location of the cool edge of the instability strip. 



the CoRoT fiel d in the galactic anti-center, such as the F- 
star HP 49933 (|Appourchaux et a l. 2008; Bcnoma r et al.1 
120091 : iKallinger et al.l l2010a| ), which has already bee n 
subject to interferometric follow-up (Bi got et al.l 1201 ID . 
However, due to poor weather conditions during the win- 
ter seasons on Mt. Wilson, not enough data has yet been 
collected for these targets. 

In the remainder of this section we summarize the ba- 
sic parameters of our target stars derived using classical 
methods and measurements available in the literature. 
Table Q] lists the complete target sample of our study, 
with spectral types taken from the HP catalog. Nine 
of the ten stars in our sample have atmospheric parame- 
ters derived from modeling several hundred lines in high- 
resolu tion spectra using t he VWA packag e 
I2010D . as presented by iBrunttj (2009), 



Bruntt et al 



. .... , ... .. .Bruntt et al., 

2012) and lThygesen et al.l ([20121 ). These are also listed 
in Table [1] For HP 189349, we have analyzed a spec- 
trum obtained with the NARVAL spectrograph at the 
Pic du Midi Observatory using three different methods : 
VWA ([Bruntt et al.ll27)Toh . ROTFI T dFrasca et al II2TM 
20061) and the met hod d escribed bylSantos et al.l ([2004) 



Sousa et al.l ([2006D and iSousa et al.1 (|2008D . In two of 



the three methods the surface gravity was fixed to the 
value calculated from asteroseismic scaling relations (see 
next section). The resulting spectroscopic parameters 
for each method are listed in Table [2j and we have 
adopted a weighted mean of all three methods, given 
in Table [TJ for the remainder of this paper. Note that 
for HP 17370 1 spectroscopic pa r amete r s have also been 
publis hed bvlValenti k, Fischerl (120051) . iMishenina et al.l 
(2004) and iKovtvukh et al.l ([20031 ). which are also listed 
in Table [2] for comparison. The published values are in 
good agreement with the values adopted here. 

All st ars in our sample h ave measured Hipparcos par- 
allaxes ([van Leeuwenl [2007[ ) . with uncertainties ranging 
from ~ 1 to 10%. All unevolved stars in our sam- 
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TABLE 1 

Fundamental properties of target stars using available literature information. Stars are separated into main-sequence 

and subgiant stars (top) and red giants (bottom). 



HD 


KIC 


Sp.T. 


y 


B — V 


Spectroscopy 


Hipparcos 












-L eft 


l°g 9 


[Pp /HI 






173701 1 


80061 61 


KOV 


7.514 


0.878 


5390(60) 


4.49(3) 


+0.34(6) 


37.47(49) 


0.64(2) 


175726 2 




G5V 


6.711 


0.571 


6070(45) 


4.53(4) 


-0.07(3) 


37.73(51) 


1.17(4) 


177153 1 


6106415 


GOV 


7.205 


0.558 


5990(60) 


4.31(3) 


-0.09(6) 


24.11(44) 


1.84(7) 


181420 2 




F2V 


6.561 


0.434 


6580(105) 


4.26(8) 


+0.00(6) 


21.05(48) 


4.1(2) 


182736 1 


8751420 


GOIV 


7.022 


0.800 


5264(60) 


3.70(3) 


-0.15(6) 


17.35(41) 


4.9(2) 


187637 1 


6225718 


F5V 


7.520 


0.500 


6230(60) 


4.32(3) 


-0.17(6) 


19.03(46) 


2.1(1) 


175955 3 


10323222 


KOIII 


7.014 


1.171 


4706(80) 


2.60(1) 


+0.06(15) 


7.62(38) 


38(4) 


177151 3 


10716853 


KOIII 


7.040 


0.994 


4898(80) 


2.62(1) 


-0.10(15) 


4.92(38) 


72(12) 


181827 3 


8813946 


KOIII 


7.188 


1.012 


4940(80) 


2.81(1) 


+0.14(15) 


4.23(43) 


83(18) 


189349 


5737655 


G5III 


7.305 


0.878 


5118(90) 


2.4(1) 


-0.56(16) 


5.32(47) 


49(9) 



B and V magnitudes are Tycho photometry (IPerrvma n & ESA 1997]) converted into the Johnson system using the calibration 

by Bcsscll (2000). Spectroscopic parameters were adopted from ^Bruntt et al.l (|2012l ). a jBruntt (200!|) and J^fhygcscn et al.l 
(2012). Spectroscopic parameters for HD 189349 are the weighted average of three results presented in this work (see text and 
Table 0. Brackets indicate the uncertainties on a parameter (note that this notation has been adopted throughout the paper). 



pie are at distances < 60 pc and hence reddening is 
expected to be negligibl e (see iMolenda-Zakowicz et al.l 
l21J09t IBruntt et alj |2012f >. Hence, we assumed zero red- 
dening for all unevolved stars with an uncertainty of 
0.005 mag. For the giants, we have estimated redden- 
ing by comparing observed colors to synthetic photome- 
try of models matching the spectroscopic parameters in 
Table [TJ as described in more detail in Section 13.31 To 
estimate an uncertainty, we have compared these values 
to E( B — V) values listed in the Kepler Input Catalog 
(KIC, 1 rown et aLl 120111) for nearby st ars and to esti- 
mates from the 3-D extinction model by iDrimmcl et al.l 
(2003). The mean scatter between these methods for all 
stars is 0.02 mag, which we adopt as our uncertainty in 
E(B — V) for the giants in our sample. Finally, we used 
the spectroscopically determined effective temperatures 
and metallicities to estimate a bolo metric correctio n for 
each star using the calibrations by iFlowerl (|1996l) and 
lAlonso et al.l (11999|) wit h appropriate zero-points as dis- 
cussed in lTorres n 2010D . yielding the stellar luminosity 
given in the last column of Table [1] Figure 1 shows an 
H-R diagram of our target stars, according to the prop- 
erties listed in Table [TJ to gether with solar-meta llicity 
BaSTI evolutionary tracks ([Piet rinfcrni et al. 2004]). 

3. OBSERVATIONS 

3.1. Asteroseismology 

The asteroseismic results presented in this paper are 
based on observations obtained by the Kepler and 
CoRoT space telescopes. Both satellites deliver near- 
uninterrupted, high S/N time series which are ideally 
suited for asteroseismic studies. In this paper, we fo- 
cus on two global parameters: the frequency of max- 
imum power (i-'max) and the large frequency separa- 
tion (A^). These are frequently used to determine 
fundamental proper ties of main-sequence and red gi- 
ant st ars (see, e.g.. iStello et al.l 12009b!: iKallinger et all 
MObTlcl: IChaplin et al.1 120111: iHekker et all I20il11rl 
Silva Aguirre et al.ll201lt Irluber et al.ll2011l) . For a gen- 
eral introduction to solar-like oscillations, we refer the 



TABLE 2 

Atmospheric parameters for stars from different methods. 



HD 


KIC 


T c ff 


logff 


[Fe/H] 


Ref 


173701 


8006161 


5390(60) 


4.49(3) 


0.34(6) 


1* 






5399(44) 


4.53(6) 


0.32(3) 


2 






5423(20) 


4.4(2) 


0.2(1) 


3 






5423(10) 






4 


189349 


5737655 


5070(100) 


2.4(1) 


-0.7(1) 


5a* 






5145(63) 


2.4(1) 


-0.54(5) 


5b* 






5163(71) 


2.9(2) 


-0.44(11) 


5c 



* lo g g fixed to asteros ei smic value. (1) Bruntt c t al.1 ( 20121 ) . 
(2) IValenti fc Fischer! (12001 ). (3) IMishenina eTaLI (|200l ). 
4) iKovtvukh et al l l|2003l), (5) this paper : a - VWA 
Bruntt et alj 120101), b - ISantos et al.1 (|2004l ) ISousa et al 
20061 ) and ISousa et al.1 (|2008l ). c - ROTFIT (|Frasca et al 



120031 12006). 



reader to the review bv lBeddind ff20lTh . 

Figure [5] presents the power spectrum for each star, 
sorted by the frequency of maximum power (v m ax)- In 
most cases, a clear power excess due to solar-like oscil- 
lations is visible. A summary of the datasets used in 
our analysis, as well as the derived asteroseismic param- 
eters, is given in Table [3J The a nalysis of Kepler stars 
is based on eithe r short-cadence dGilliland et aLll2010al) 
or long-cadence ([Jenkins et al.ll2010D data up to Q10, 
which were corre cted f or instrumental trends as described 
(l20Tll) . 



in iGarria et al.l (|2011l ). Global asteroseismic parameters 
were extract e d usin g the automated analysis pipeline by 
iHuber et al.1 ([2009D , which has been shown to agree well 
with other methods (|Hekker et al.l l2011ct iVerner et al.l 
120111 ) - Due to the length and very high S/N of the Kepler 
data, the modes are resolved and uncertainties on z^ max 
and (particularly) are dominated by the adopted 
method (e.g., the range over which /S.v is determined) 
rather than measurement errors. To account for this, we 
added in quadrature to the formal uncertainties an uncer- 
tainty based on the scatte r of different methods used by 
iSilva Aguirre et ahl (|2012f ) for short-cadence data and by 
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Fig. 2. — Power density spectra for all stars in our sample, sorted by the frequency of maximum power (v m ax). Note the change in x-axis 
scale for main-sequence (left column), subgiant (top right column) and red-gian t stars (four bottom r ight panels). Note that the high peak 
at ~4500/^Hz for HD 187637 is a known artefact of Kepler short-cadence data (Gilliland et al. 2010a). 



IHuber et al.l (120111) fo r long-cadence data. The analysis 
by IHuber et al.l (|2011[ ) was based on data spanning from 
QO-6, which in most cases was sufficient to resolve the 
mode s and reliably estimate f ma x and Av (|Hekker et al.l 
2012). In general, the uncertainties on the asteroseismic 
parameters for most Kepler stars are negligible compared 
to the uncertainties on other observables. A notable ex- 
ception is HD 189349, with a relatively large uncertainty 
of ~ 4% in the large frequency separation. Inspection of 
the power spectrum shows that the modes for this star 
are very broad, making a determination of Av difficult. 
We speculate that the unusually broad modes may be 
related to the low metallicity of this object, but a more 



in-depth analysis is beyond the scope of this paper. 

For the two CoRoT stars in our sample, we have re- 
analyzed publicly avai l able data using the method de- 
scribed in lHuber et all (|2009| ). Our results for HD181420 
are in good a greem ent with the values published by 
iBarban et ail ([20091) . For HD175726, our analysis did 
not yield significant evidence for regularly spaced peaks, 
and yielded only marginal evidence f or a power excess 
at 1900 ± 200^iHz. iMosser et all (|2009f ) have argued that 
this power excess is compatible with solar-like oscillations 
and showed evidence for a large variation of Av with fre- 
quency, which could be responsible for the null-detection 
in our analysis. We have adopted the published value 
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TABLE 3 

ASTEROSEISMIC OBSERVATIONS AND MEASURED PARAMETERS. 



11 L) 






i (a) 


uuty cycle \/o) 






I/O ( Ul 


onnc i/?i 


tv",,,-- l ori r\c in 

ivepler t^o-lU 


KC'7 
00 I 


no 
Uz 


ooly(yo J 


1 A(\ Q ( A \ 

14y.o(4j 


HOt ZD 




OOIXO _L QJTtCUl 


Zo 


on 


1 Q1 ^900^ 
lyio^zuu ) 


Q7 

y i .z^o ) 


177153 


6106415 


Kepler SC Q6-8.10 


461 


70 


2233(60) 


104.3(3) 


181420 




CoRoT LRcOl 


156 


90 


1574(31) 


75.1(3) 


182736 


8751420 


Kepler SC Q5,7-10 


557 


75 


568(15) 


34.6(1) 


187637 


6225718 


Kepler SC Q6-10 


461 


91 


2352(66) 


105.8(3) 


175955 


10323222 


Kepler LC Q0-10 


880 


91 


46.7(1.1) 


4.86(3) 


177151 


10716853 


Kepler LC Q 1-7,9- 10 


869 


77 


48.8(1.1) 


4.98(7) 


181827 


8813946 


Kepler LC Ql-10 


869 


91 


73.1(1.2) 


6.45(7) 


189349 


5737655 


Kepler LC Ql-10 


869 


91 


29.9(1.1) 


4.22(16) 



1 Detection adopted from Mosser et al. (2009). 



for Av by iMosser et al.l (|2009f ) and a value for v max cor- 
responding to the maximum of the power excess in the 
spectrum, with a conservative uncertainty of 10%. 

3.2. Interferometry 

Interferometric observations were made with the Pre- 
cision Ast ronomical Visible O bservations (PAVO) beam 
combiner (Irel and et alj|2008l) at the Center for High An- 
gular Resolution Astrono my (CHARA) on Mt. Wilso n 
observatory, California (jten Brummelaar et al.l [2005). 
Operating at a central wavelength of A = 0.7 /im with 
baselines up to 330 m, PAVO at CHARA is one of the 
highest angular-resolution instruments world-wide. 

A complete d e script ion of the instrument was given 
bv llreland et al . (2008), and we summarize the basic as- 
pects here. The light from up to three telescopes passes 
through vacuum tubes and into a series of optics to com- 
pensate the path difference. The beams are then colli- 
mated and passed through a non-redundant mask which 
acts as a bandpass filter, and spatially modulated in- 
terference fringes are formed behind the mask. The in- 
terference pattern is then passed through a lenslet ar- 
ray and a prism, producing fringes in 16 segments on 
the CCD detector, each being spectrally dispersed in 
several independent wavelength channels. Major advan- 
tages of the PAVO design are high sensitivity (with a 
limiting magnitude of R ~ 8 mag in typical seeing con- 
ditions), increased information through spectral disper- 
sion, and high spatial resolution through operating at vis- 
ible wavelengt hs. First PAVO sc ie nce results have been 
pres ented by iBazot et al.l (|201lD . iDerekas et all (|201lD 
and lHuber et alj (|2012f T" For our analysis, we have used 
PAVO observations in two-telescope mode, with base- 
lines ranging from ~ 1 10 — 330 m. 

Interferometric observations require careful calibration 
of the observed visibilities. Ideally, this is achieved by 
observing bright, unresolved point sources as closely as 
possible to the target object in time and distance. For 
PAVO observations of targets as small as in our case, this 
means calibrating with late B to early A stars since at the 
PAVO magnitude limit these stars are distant enough to 
have significantly smaller diameters (0.1-0.15 mas) than 
our target stars. Table Q] lists all calibrators that were 
used in our analysis. Expected sizes are ca lculated using 
the V — K relation bv lKervella et al.l (|2004f ) for dwarf and 
subgiant stars. V band magnitudes have been taken from 



TABLE 4 

Calibrators used for interferometric observations. 



HD 


Sp.T. 


V-K 


E(B - V) 


0v-K 


ID 


171654 


A0V 


-0.067 


0.036 


0.141 


c 


174177 


A0V 


0.249 


0.020 


0.191 


gh 


176131 


A2V 


0.345 


0.012 


0.155 


ac 


176626 


A2V 


0.084 


0.026 


0.146 


ac 


177959 


A3V 


0.451 


0.029 


0.152 


b 


178190 


A2V 


0.381 


0.027 


0.157 


bd 


179095 


A0V 


-0.069 


0.022 


0.129 


gh 


179124 


B9V 


0.280 


0.095 


0.146 


el 


179483 


A2V 


0.316 


0.028 


0.144 


e 


179733 


A0V 


0.211 


0.038 


0.117 


ac 


180138 


A0V 


0.075 


0.045 


0.128 


c 


180501 


A0V 


0.147 


0.027 


0.117 


gh 


180681 


A0V 


0.112 


0.031 


0.111 


acei 


183142 


B8V 


-0.462 


0.060 


0.093 


ei 


184147 


A0V 


0.007 


0.019 


0.121 


egi 


184787 


A0V 


0.034 


0.017 


0.154 


ef 


188252 


B2III 


-0.461 


0.047 


0.155 


ce 


188461 


B3V 


-0.461 


0.109 


0.095 


ef] 


189845 


A0V 


0.136 


0.053 


0.127 


fj 


190025 


B5V 


-0.230 


0.157 


0.084 


j 


190112 


A0V 


0.067 


0.027 


0.113 


f 



List of dropped calibrators: HD179395, HD181939, 
HD182487, HD184875, HD189253; "ID" refers to the ID of 
the target star for which the calibrator has been used (see 
column 3 of Table O . 



the Tycho catalog and were co nverted into t he Johnson 
system using the calibration by|Bessell (2000). K. magni- 
tudes were adopted from 2MASS ()Skrutskie et alT 2006). 
Interstellar reddening for eac h calibrator was estim ated 
using the extinction model bv lDrimmel et all (|2003l ). 

Although we have checked each calibrator in the lit- 
erature for possible multiplicity, rotation and variabil- 
ity prior to observations, our data show that roughly 
1/4 of all observed calibrators are more resolved than 
expected, and therefore potentially unsuitable for cali- 
bration. These calibrators are listed at the bottom of 
Tabled] Possible reasons for this include previously un- 
detected binary systems and rapid rotation causing de- 
viations from spherical symmetry. 
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TABLE 5 

Interferomethic observations and measured parameters. 



un 
riu 






Scans /Nights 


Baselines 


fJ-R 


PUD 




n 

C(V-iT) 


#IRFM 


1 ( o ( Ul 


oUUOlOl 


a 


O I A 


S2E2,S1W1,S1E2,S1E1 0.59(4) 


O Q1 A (A \ 


\J.ooZ[o ) 


\j.666\p ) 




1 ( ( zo 






1 /9 
O j Z 


S1W2,S1E1 


0.51(5) 


U.ool^O ) 




U.OOD^O ) 




177153 


6106415 


c 


7/3 


S2E2,S1W1,S1E2 


0.51(5) 


0.276(6) 


0.289(6) 


0.285(4) 


0.28(1) 


181420 




d 


5/2 


S1W2,S1E1 


0.48(5) 


0.32(1) 


0.34(1) 


0.312(5) 


0.31(1) 


182736 


8751420 


c 


5/4 


S2W2,S2E2,S1W1 


0.59(4) 


0.412(3) 


0.436(5) 


0.429(6) 


0.44(2) 


187637 


6225718 


f 


6/3 


S2E2,S1E2,S1E1 


0.49(5) 


0.222(5) 


0.231(6) 


0.222(3) 


0.22(1) 


175955 


10323222 


g 


4/2 


W1W2,S2W2 


0.67(2) 


0.634(9) 


0.68(1) 


0.70(1) 


0.66(3) 


177151 


10716853 


h 


4/2 


W1W2,S2W2 


0.64(3) 


0.541(8) 


0.57(1) 


0.57(1) 


0.53(2) 


181827 


8813946 


i 


3/2 


S2W2,S1W1 


0.64(3) 


0.443(3) 


0.473(5) 


0.516(9) 


0.49(2) 


189349 


5737655 


j 


4/3 


S2W2,S1W2 


0.58(4) 


0.399(4) 


0.420(6) 


0.444(9) 


0.44(2) 



"ID" can be used in Table [4] to identify which stars have been used to calibrate this target. All angular diameters are given in 
units of milli-arcseconds. Baselines are sorted from shortest to longest length for a given target. 



Figure |3] presents the calibrated squared-visibility mea- 
surements as a function of spatial frequency for all targets 
in our sample, with a summary of observations given in 
Table [5] We have collected at least three independent 
scans for each target over at least two different nights, 
and the visibilities of each target were calibrated with at 
least two different calibrators (see also Tabled]). Note 
that each scan typically produces a measurement of vis- 
ibility in 20 independent wavelength channels, resulting 
in a total of ^1000 visibility measurements in our cam- 
paign. 

For each target we fitted the fo llowing limb-darkened 
disc model to the observations (|Hanburv Brown et al.1 
[1971 : 



V = 



with 



MA 



I ' 3 , 

/, /lW . , /0 U/2 J 3/2(z) 

X x i ' z 



x = ttBOldX 



(1) 



(2) 

Here, V is the visibility, fix is the linear limb-darkening 
coefficient, J n (x) is the n-th order Bessel function, B 
is the projected baseline, #ld is the angular diame- 
ter after correction for limb-darkening, and A is the 
wavelength at which the observation was made. Linear 
limb-darkening coefficients in the R band for our tar- 
gets were estimated by i nterpolating the model grid of 
iClaret fc Bloem cn (201 ll) to the spectroscopic estimates 
of T e g, log g and [Fe/H] (Table [T]) for a microturbulcnt 
velocity of 2 km s~ 1 . Uncertainties on the limb-darkening 
coefficients were estim ated from the diff e rence in the 
methods presented by IClaret fc Bloemenl (|201lD . The 
choice of the limb-darkening model has little effect on 
the final fitted angu lar diameters. De t ailed 3-D hydro- 
dynam ical models bvlBigot et al.1 (120061) , IChiavassa et al.1 
(2010) and IChi avassa et al. (2012]) for dwarfs and giants 
have shown that the differences to simple linear limb 
darkening models are 1% or less in angular diameter for 
stars with near solar-metallicity. For a moderately re- 
solved star with V 2 ~ 0.5, a 1% change in angular di- 
ameter would arise from a change of less than 1% in V 2 , 



which is less than our typical measurement uncertainties. 

The procedure used to fit the model and estimate the 
uncertaint y in the derived ang ular diameters was de- 
scribed by iDerekas et al.l (|201lD . In summary, Monte- 
Carlo simulations were performed which took into ac- 
count uncertainties in the adopted wavelength calibra- 
tion (0.5%), calibrator sizes (5%), limb-darkening coef- 
ficients (see Table [5]) , as well as potential correlations 
across wavelength channels. The resulting fitted angular 
diameters of each target, corrected for limb-darkening, 
are given in Table [S] We also give the uniform-disc di- 
ameters in Table [SJ which were derived by setting fix = 
in Equation (1). 

A few comments on our derived diameters are neces- 
sary. Firstly, one calibrator in our sample (HD 179124), 
which is the main calibrator for HD 181420, was recently 
found to be a rapidly rota ting B star with vsmi = 
290km/s (jLefever et al.H2010D . This introduces an extra 
uncertainty on the estimated calibrator diameter. We 
have accounted for this by assuming a 20% uncertainty 
in the calibrator diameter, which roughly corresponds to 
the maximum ch ange in the average diameter exp ected 
for rapid rotators (jDomiciano de Souza et al.l l2002). Sec- 
ondly, a few of our target stars (e.g. HD 187637) are only 
about 50% bigger in angular size than their calibrators. 
This means that the uncertainties on the derived diam- 
eters will be strongly influenced by the assumed uncer- 
tainties of the calibrator diameters, which in our case 
are 5%. While such an uncertainty is reasonable com- 
pare d to the scatter in th e photometric calibrations (see, 
e.g.. iKervella et al~ll2004h . the diameter measurement it- 
self will only be scientifically useful if the uncertainty 
in the measured diameter is smaller than the precision 
of indirect techniques. Further data at longer baselines 
with smaller calibrators will be needed to reduce the un- 
certainties for these targets. 

Indirect techniques to estimate angular dia meters in- 
clude surface brightn e ss rel ations (see, e.g., Ivan Belle] 
119991: IKervella et all 120041) and the infrared flux 
method (IRFM, see, e .g., iRarmrez fc Melendezl 120051: 
ICasagrande et al.l 120101 ) . Figure H compares our mea- 
sured angular diameters with predictions using the 
(V — K) surface-brightn e ss rela tion for dwarfs and sub- 
giants by IKervella et~aTI (|2004D and the IRFM method 
coupled with asteroseism ic constraints, as described in 
iSilva Aguirre et al.l (|2012f ). For the (V - K) relation we 
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Fig. 3. — Squared visibility versus spatial frequency for all stars in our sample. Red solid lines show the fitted limb-darkened disc model. 
The order of panels is the same as in Figure [2] Note that the error bars for each star have been scaled so that the reduced x 2 equals unity. 



have adopted a 1% diameter uncertainty for all stars 
(jKervella et al.ll2004r i. We find good agreement for all 
stars for both methods, with a residual mean of —2 ±2% 
and +2 ± 2% for (V - K) and IRFM, respectively, both 
with a scatter of 5%. Our results therefor e seem to con- 
firm that the relation by iKervella et al.1 (120041) is als o 
valid for red giants, as sugg ested bv iPiau et al l (|2011h . 
and that combining the IRFM method with astero seismic 
constraints, as done bv lSilva Aguirre et al.l ()2012[ ). yields 
accurate diameters for both evolved and unevolved stars. 

These tests of indirect methods are encouraging. We 
emphasize that interferometry remains an important tool 
to validate these methods for a wider range of evolution- 
ary states, chemical compositions, and distances. The 



(V — K) relation, for example, is based on an empirical 
relation calibrated using nearby stars that does not take 
into account potential spread due to different chemical 
compositions, and is only valid for de-reddened magni- 
tudes. An illustration of the importance of using interfer- 
ometry is HD 181827, which shows a significantly smaller 
measured diameter than predicted from (V — K). This 
smaller diameter is also in agreement with asteroseismic 
results, which suggest a smaller radius (see Section 4.1). 

3.3. Bolometric Fluxes 

To estimate bolometric fluxes for our target sample, 
we first extracted synthetic fluxes fr om the MARCS 
database of stellar model atmospheres ()Gustafsson et al.l 
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Fig. 4. — Fractional differences between angular diameters mea- 
sured with PAVO to diamet ers determined using the (V — K) 
surface-brightness relation by Kcrvclla ct al. ( 2004) (upper panel) 
and using the i nfrared flux method with asteroseismic constraints, 
as described in Silva Aeuirrc ct al. (2012) (lower panel). Black di- 
amonds show main-sequence and subgiants stars, and red triangles 
show giant stars. HD numbers of each target are labelled in the 
upper panel. 



2008). We used models with standard chemical com- 
position, with the microturbulence parameter set to 
lkms -1 for plane-parallel models (unevolved stars) 
and 2kms~ 1 for spherical models with a mass of 
IMq for red giants. We then multiplied the syn- 
thetic stellar fluxes by the filter responses for the 
Johnson-Glass-Cousins UBVRIJHKL, Tycho B T V T 
and 2MASS JHK S systems and integrated the result- 
ing fluxes to calculate synthetic magnitudes for each 
MARCS m odel. Filter responses a nd zeropoints were 
taken fromlBessell fe Murphvl (120121) (UBVRI, ByVr) , 
ICohen et all (120031 ) (2MASS1. and iBessell et all pi ) 
(JHKL). We note that synthetic photometry calculated 
using MARCS models has previously bee n validated us- 
ing observed colors in stellar clusters dBrasseur et al.l 
[20101: iVandenBerg et al.|[20Tot) . To check the influence 
of the chosen mass for the spherical models, we have re- 
peated the above calculations for typical red giant models 
with T cff = 5000 K and logg = 2 - 3. The fractional dif- 
ferences in the integrated flux for each filter for masses 
ranging from 0.5 — 5M Q was found to be less than 0.5% 
in all bands, and are therefore negligible for our analysis. 

The amount of photometry in the literature for our 
sample is unfortunately small. The targets are generally 
too faint to have reliable magnitudes in the Johnson- 
Glass-Cousins system, and they are too bright to have a 
full set of SDSS photometry in the KIC. To ensure consis- 
tency of our bolomctric fluxes we only used photometry 
that is available for all stars in our sample, namely Ty- 
cho2 B T V T and 2MASS JHK S magnitudes. The adopted 
photometry and uncertainties are listed in Table [5] 

To calculate bolom etric fluxes, we large ly followed the 
approach described in lAlonso et al.l |l995). For each tar- 
get star, we first found the six models bracketing the 
spectroscopic determinations T c g, \ogg and [Fe/H], as 




500 



1000 1500 2000 

Wavelength (rim) 



i50(; 



Fig. 5. — Spectral energy distibutions of Procyon, the Sun and 
Arcturus to test our method to determine bolometric fluxes. Black 
lines are MARCS models with parameters as given in each panel, 
smoothed with a constant spectral resolution A/ A A ~ 200 (cor- 
responding to a width of ~2.5 nm in the V-band). Observed and 
synthetic BVJHK photometry is shown as red diamonds and blue 
squares, respectively. All fluxes have been normalized to 1 in the 
band with highest flux for a given star. The inset shows the re- 
sult of Monte-Carlo simulations to estimate the uncertainty in the 
bolometric flux as described in the text. Dashed and dotted lines 
show the literature values and l-cr uncertainties. 



given in Table [TJ We then transformed the synthetic 
BtVtJHK s magnitudes of each model into fluxes, and 
numerically integrated these fluxes using the pivot wave- 
len gth for each filter respon se, calculated as described 
by iBessell fc Murphvl ([2012) (note that this choice of a 
reference wavelength is independent of the spectral type 
considered). The numerical integration yielded an esti- 
mate /i n t, which we then compared to the true bolometric 
flux, /boi = cr^eff) where a is the Stefan-Boltzmann con- 
stant. This yielded a correction factor c = /i n t//boi f° r 
each of the six models, which gave the percentage of flux 
included when integrating the photometry over discrete 
wavelengths. The final bolometric flux was then calcu- 
lated by integrating the observed fluxes the same way 
as the model fluxes, and dividing the resulting estimate 
by the correction factor c found by interpolating the six 
corrections factors to the spectroscopic estimates of T Q g , 
logg and [Fe/H]. Note that this interpolation was neces- 
sary because the step size of the MARCS grid is typically 
larger than the uncertainties of the spectroscopic param- 
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Fig. 6. — Spectral energy distributions for all stars in our sample. Black lines are the MARCS models with parameters as given in each 
panel, smoothed with a constant spectral resolution A/AA ~ 200 (corresponding to a width of ~2.5nm in the V-band). Red diamonds and 
blue squares show the observed and model fluxes in the BtVtJHK bands, respectively. All fluxes have been normalized to 1 in the band 
with highest flux for a given star. The order of panels is the same as in Figure [2] 



eters. Uncertainties in the derived bolometric fluxes were 
found by perturbing the input photometry and the spec- 
troscopic parameters according to their estimated uncer- 
tainties (see Tables Q] and [6]), repeating the procedure 
5000 times, and taking the standard deviation of the re- 
sulting distribution. 

To test this approach, we have used the same method 
for three bright stars that span a similar range of evolu- 
tionary stages as our sample and for which bolometric 
fluxes have been well determined: Procyon, the Sun, 
and Arcturus. Since Tycho and 2MASS photometry 
is not available for such bright stars, we have used 
BVJHK photometry to mimic the available informa- 



tion for our target sample. Photometry has been taken 
from the General Cata log of Photometric Data (GCDP, 
iMermilliod et all I1997D f or Procyon and Arcturus, and 
from iColina et al .1 (j!99ffl for the Sun. Figure [5] shows 
the spectral energy distributions (SEDs) of all three 
stars, comparing the MARCS model that best matches 
the physical parameters of each star (black line) to the 
observed and synthetic fluxes in the BVJHK bands 
(red and blue squares, respectively). Note that the 
MARCS models have been smoothed to a spectral res- 
olution of A/AA ~ 200 for better visibility. The insets 
show the distributions of the Monte-Carlo simulations 
described above compared to the literature values of the 
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TABLE 6 

Broadband photometry, estimated reddening and bolometrio fluxes for all target stars. 



HD 


KIC 




V T 


J 


H 


K 


E(B - V) 


Fbol 


(10~ s erg s~ 


1 cm ' 2 ) 


















MARCS 


ATLAS+9 


ATLAS+Vt 


173701 


8006161 


8.606(16) 


7.610(11) 


6.088(21) 


5.751(16) 


5.670(21) 


0.000(5) 


2.89(5) 


2.93(4) 


2.86(3) 


175726 




7.401(15) 


6.780(10) 


5.703(24) 


5.418(34) 


5.346(20) 


0.000(5) 


5.4(1) 


5.40(7) 


5.41(6) 


177153 


6106415 


7.872(15) 


7.275(11) 


6.145(20) 


5.923(26) 


5.829(23) 


0.000(5) 


3.39(7) 






181420 




7.059(15) 


6.604(10) 


5.748(21) 


5.560(33) 


5.513(26) 


0.000(5) 


6.0(2) 


6.2(1) 


5.99(7) 


182736 


8751420 


8.021(16) 


7.103(10) 


5.515(24) 


5.135(27) 


5.028(16) 


0.000(5) 


4.77(8) 


4.74(5) 


4.67(5) 


187637 


6225718 


8.118(15) 


7.580(11) 


6.544(21) 


6.346(29) 


6.283(18) 


0.000(5) 


2.55(5) 






175955 


10323222 


8.513(16) 


7.146(10) 


4.999(24) 


4.442(31) 


4.318(17) 


0.09(2) 


7.2(3) 






177151 


10716853 


8.281(16) 


7.144(10) 


5.264(18) 


4.820(33) 


4.686(20) 


0.04(2) 


5.6(3) 






181827 


8813946 


8.476(15) 


7.300(10) 


5.453(35) 


4.997(15) 


4.872(21) 


0.04(2) 


4.8(2) 






189349 


5737655 


8.411(16) 


7.411(10) 


5.638(24) 


5.181(21) 


5.124(29) 


0.07(2) 


4.6(2) 







Tycho2 and 2MASS photometry are taken from IHgg eFaLI ||2000T) and lCutri eTaT] {2003). 



bolometric flux (dashed line) and their 1-a uncertain- 
ties (dotted lines ). Literature bolometric fl uxes have 
been tak en from iRamrrez fc Allende Prietol (1201 ID for 
Arcturus, [Aufdenb erg et al.l (|2005l ) and iFuhrman n et al.l 
(fl99l for Procyon, and we have adopted an effective 
temperature of 5777±3K for the Sun. In all three cases, 
the bolometric flux using our method is recovered within 
1-a, with a maximum deviation of ~0.5 a for Arcturus. 

Figure [5] shows the SEDs of all target stars with the ap- 
propriate models for each star, and Table[6]lists our bolo- 
metric fluxes based on the procedure described above. 
We note that for the red giants in our sample, inter- 
stellar reddening cannot be neglected. To estimate red- 
dening using the SED, we adopted the reddeni ng law by 
lO'Donnelll (|199l (see also lCardelli"eTaT][r989l ) and iter- 
ated over E(B — V) to find the observed colors that best 
fit the colors of the six models bracketing the spectro- 
scopic parameters in Tabled] We then again interpolated 
to the spectroscopic T e g, \ogg and [Fe/H] values, anal- 
ogously to the correction factor described above. The 
derived reddening estimates for the giants are listed in 
Table [U 

To further test these results, we have used an in- 
dependent method to determine bolometric fluxes for 
four stars by com bining publicly availa ble flux-calibrated 
ELODIE spectra (jPrugniel et al.ll2007D. broadband pho- 
tome try and ATLAS9 models ( |Castelli fc Kuruczl 120031 
I2004D . We started by calculating a grid of ATLAS9 mod- 
els in the 3-cr error box of the spectroscopically deter- 
mined T e g, logg and [Fe/H] (see Table [T]). Each model 
spectrum was then multiplied by the BtVtJHK s filter 
passbands and integrated over all wavelengths to com- 
pute a synthetic flux in each band. Model fluxes were 
then calibrated into fluxes received on Earth using ei- 
ther the measured angular diameter or the Tycho Vt 
magnitude. To find the model that best fits the photo- 
metric data we then compared the grid of model fluxes 
with the observed fluxes, calculated using the same ze- 
ropoints as in the procedure described above. Finally, 
the bolometric flux of each star was determined by inte- 
grating the ELODIE spectrum between 390 and 680 nm 
together with the synthetic ATLAS9 model (covering the 
wavelength ranges < 390 nm and > 680 nm) that best fits 
the observed photometry 

To estimate uncertainties the above procedure was re- 
peated 100 times, drawing random values for the ob- 



served photometry given in Table[6j and adding the stan- 
dard deviation of the resulting distribution in quadra- 
ture to the uncertainty of the total flux of the ELODIE 
spectra. The final values for the two different cali- 
bration methods are given in Table [6] The derived 
bolometric fluxes agree well with the estimates from 
MARCS models, reassuring us that the model depen- 
dency and adopted method have little influence com- 
pared to the estimated uncertainties. We note that we 
have also compared our bolometric fluxes with estimates 
derived from the infrared flux method, as described in 
iSilva Aguirre et al.1 (|2012l ). A gain, we have found good 
agreement with our estimates within the quoted uncer- 
tainties. 

4. FUNDAMENTAL STELLAR PROPERTIES 

4.1. Asteroseismic Scaling Relations 

The large frequency separation of oscillation modes 
with the same spherical degree and consecutive ra dial or- 
der is closely related to the mean density of star (|Ulrichl 
[l986T k 

Au oc M 1/2 iT 3/2 . (3) 

Additionally, iBrown et al.l (|1991|) argued that the fre- 
quency of maximum power (f ma x) for solar-like stars 
should scale with the acoustic cut-off f requency, which 
was used bv iKjeldsen k. Bedding! (|1995f ) to formulate a 
second scaling relation: 

*/ max oc MR- 2 Tj /2 . (4) 

Provided the effective temperature of a star is known, 
Equations (J3j> and (jlj allow an estimate of the stellar 
mass and radius. This can be done by either combin- 
ing the two equations (the so-called direct method, see 
iKallinger et alJl2010d ) or by comparing the observed val- 
ues of Ai/ and f max with values calculated from a grid 
of evolutionary models (the so-called gr i d-based method , 
see iStello et al.ll2009bt iBasu et al.|[2010l: iGai et al.|[201lT) . 

Our intcrfcromctric observations, presented in the Sec- 
tion 3.2, allow us to test Equations ([3]) and (|4]). Using the 
Hipparcos parallaxes in combination with the angular di- 
ameters, we have calculated linear radii for our sample 
of stars, which are listed in Table [7J These are compared 
to asteroseismic radii calculated using Equations ([3]) and 
((U (using T e ff values taken from Table [J) in Figure [7j 
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Fig. 7. — Comparison of stellar radii measured using interfer- 
ometry and calculated using asteroseismic scaling relations. Black 
diamonds show our Kepler and CoRoT sample, and red asterisks 
show several bright stars as indicated in the plot for comparison. 
The dashed line marks the 1:1 relation. 
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Fig. 8. — Comparison of !^ max measured from asteroseismology 
and calculated using independent measurements of R, M and T cf f . 
Black diamonds show the Kepler and CoRoT sample, and red as- 
terisks show several bright stars as indicated in the plot for com- 
parison. The dashed line marks the 1:1 relation. 



Note the influence of T e g on Equation (|4]) is small: for 
solar T e ff a variation of 100 K causes only a 0.9% change 
in !/ max , which is significantly smaller than our typical 
uncertainties (see Table [3]) . 

The comparison in Figure [7] is very encouraging, show- 
ing an agreement between the two methods within 3-cr 
in all cases. The overall scatter about the residuals is 
~13%, and we do not observe any systematic trend as 
a function of size (and therefore stellar properties). We 
note that two of the stars in our sample ( HP 173701 and 
HP 17 7153) have also been analyzed by iMathur et al.l 
(2012), who used both a grid-based approach as well as 
detailed modelling of individual oscillation frequencies to 
derive stellar radii and masses. In both cases, the radi i 
from different models presented in IMathur et al.l (|2012l ) 
slightly improve the difference to the interferometrically 
measured radius, with minimum differences of +0.4 a and 
+0.8 a compared to differences of -0.6 a and -1.0 a from 
the direct method, respectively. 

For comparison, Figure [7] also shows examples of 
bright stars for which well constrained asteroseismic 
and interferometric parameters are availabl e. We have 
adopte d values for Av and z/ max from iStello et"al"1 
(2009a) and references therein, with uncertainties 
fixed to typical values of 1% in Av and 3% in v ma , x - 
Asteroseismic observations have been o btained from 
the MOST space telesco pe for eOph (jBarban et al.l 
I2007t iKallinger etldl 120081). the CoR oT space-telescope 
for HP 49933 (jAppourchaux et all 120081) . and from 
ground-based Poppler observat i ons f o r the remain- 
ing sample (|Carrier fc Bourbanl 120031 : Kicldsc n et al.l 



pOOl IBeddine et al.1 l200l ICarrier et al.1 I2005allbl: 
Kieldsen et alJl2005i .iBedding et alJl200ft lArentoft et al.l 
120081 : iTeixeira et alJ 12009b IBazot et all 120111 ). Angular 
diameters and e ff ective tempe ratures were take n from 
iMazumdar eTaTl (120091) and iPe Ridder eFaTl (12001) 
fore Oph. IBazot et al.f (120111) for 18 Sco. iBigot et al.1 
poTfl) for HP 49933 and from lBruntt et al.l d2010D ~and 
references therein fo r the remaini n g sam ple. Parallaxes 
were adopted from Ivan Leeuwen (2007), except for a 
Cen A and B for which we have adopted the value by 
ISoderhielml (fl999l) . Figure [7] again shows agreement 
within 3-cr in all cases. Excluding HP 175726 from our 
sample due to large uncertainties in the asteroseismic 
observations, the residual scatter between asteroseismic 
and interferometric radii is 4% for dwarfs and 16% for 
giants, with mean deviations of —1 ± 1 % and +6 ± 4 %, 
respectively. This is consistent with our observational 
uncertainties and hence empirically confirms that, at 
least for main-sequence stars, asteroseismic rad ii from 
scalin g relations are accurate to < 4%. Note that[Miglio 
(2012) has previously found a similar good agreement 
for a sample of nearby stars, with a residual scatter of 
6%. 

It is well known that the scaling relation for Av 
is on more solid ground than the scaling relation 
for i/ mH , which only recently ha s been studied in 
more detail obseryationally (see, e.g.. | Stello et al.l f2009a: 



IMosser et ail 120101: IWhite et al.l |2011J) and theoretically 
(|Belkacem et al.l 1201 ID . To test Equation @, we can 



combine Equation (J3]) with the interferometrically mea- 
sured radii to calculate stellar masses, and combine these 
with T e ff to calculate f max - We compare these with the 
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measured values in Figure |S] We again observe good 
agreement within the error bars, with no systematic de- 
viation as a function of evolutionary status. Figure[8]also 
displays a comparison with measured values for a sam- 
ple of bright stars, again showing good agreement with 
our results fo r the K epler and CoRoT sample. We note 
that iBeddinel (|2011l) has shown a similar comparison for 
bright stars, and noted a potential breakdown of the ^ max 
relation for low-mass stars with f max > 4500/iHz. Since 
none of the stars in our sample have f max > 4000/iHz, 
we are unable to test this claim in our study. 

The large error bars for some stars in Figures [7J and 
[5] may cast some doubt about the usefulness of inter- 
ferometry to test scaling relations. Indeed, for the red 
giants in our sample the uncertainty in the interferomet- 
ric radius is completely dominated by the uncertainty 
in the parallax. For these stars the PAVO data will be 
most valuable to measure the effective temperature by 
combining the angular diameter with an estimate of the 
bolometric flux, which can then be compared to indi- 
rect T e ff estimates from broadband photometry and spec- 
troscopy (see next section). For most unevolved stars 
in the Kepler /CoRoT sample, our current uncertainties 
in the angular diameters arc comparable to the paral- 
lax uncertainties. The bright star comparison sample, 
on the other hand, is dominated by the uncertainties in 
the asteroseismic observables, which are much more diffi- 
cult to constrain from the ground or using smaller space 
telescopes. The fact that the asteroseismic uncertainties 
are almost negligible for the Kepler /CoRoT sample ex- 
plains the somewhat counter-intuitive observation that 
the error bars in Figures [7J and [5] are similar for some 
stars of the Kepler sample and for stars which are up to 
8 magnitudes brighter. This comparison underlines the 
importance of obtaining precise asteroseismic constraints 
on bright stars for which constraints are available from 
independent observational techniques. 

4.2. Spectroscopic and Photometric Temperatures 

The measurement of the angular diameter #ld of a 
star combined with an estimate of its bolometric flux /boi 
allows a direct measurement of the effective temperature: 

where a is the Stefan-Boltzmann constant. We have used 
our measured angular diameters presented in Section 3.2 
together with the bolometric flux estimates presented in 
Section 3.3 to calculate effective temperatures for our 
sample, which are listed in Table [Jj 

The model dependency of effective temperatures cal- 
culated using Equation ([5]) is small, and hence such 
estimates are important for calibrating indirect pho- 
tometric estimates such as the infrared flux method 
(see, e.g., iCasagrande et al.ll2010D . as well as spectro- 
scopic determinations for which usually strong degen- 
eracies between T K ff, logg and [Fe/H] exist (see, e.g., 
iTorres et all |2012[ ). Figure M compares the measured 
effective temperatures in our sample to estimates from 
high-resolu tion spectroscopy ( mostly using the VWA 
package bv lBruntt et al.l (|2010h . see Tabled! as well as 
photo m etric calibrations taken from ICasagrande et ah 
i|2010t) . iRami'rez fc Melendezl (|2005| ) and IBruntt et al. 




a, : 
^ -S00 7 

3 -400 r 




4500 5000 5500 6000 6500 7000 

^eff.PAVO (K) 

Fig. 9. — Comparison of effective temperatures derived in this 
study with spectroscopic estimates (panel a) and several photo- 
metric calibrations based on V — K B (panels b-d). Black diamonds 
are main-sequence sta rs, while red triang l es sho w red giant s. Note 
that t he relations by Casagrandc et al. ( 2 0101) and [Bruntt et al. 
(2012) are calibrated for main-sequence stars. 

(2012). We have chosen V — K s to calculate photomet- 
ric temperatures since this index usually gives the lowest 
resid uals as a temperature in dicator for cool stars (see, 
e.g.. ICasagrande et allfeOlCj) . The comparison in Fig- 
ure [5k< shows good agreement of our temperatures with 
spectroscopy, with a residual mean of —48 ± 39 K with 
a scatter of 124 K for all stars, and —22 ± 33 K with a 
scatter of 97 K when excluding the F-star HD 181420 for 
which the angular diameter is not well determined. We 
note that this agreement is only slightly worse (with an 
inc reased scatter by ab out 10 K) if we use the T e g va lues 
by IBruntt et all (|2012fl and iThygesen efail (|2012D for 
which no asterosei smic constraints on \ogg were used. 
IBruntt et al.l (|2010f ) noted that a slight bias for spectro- 
scopic temperatures to be hotter than intcrfcrometric es- 
timates by ~ 40 K for a sample of nearby stars, which is 
somewhat confirmed by our results, although the scatter 
is significantly larger. Our result confirms that a combi- 
nation of spectroscopy and asteroseismology can be ap- 
plied for the accurate characterization of temperatures, 
radii and masses of much fainter stars, e.g. exoplanet 
host stars observed by the Kepler mission. 
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The photometric estimates shown in Figure |Hb, [St 
and [9(1 show slight systematic d eviations. The calibra- 
tion by iCasagrande et al.l (|2010f) shows the best agree- 
ment, with only the coolest red giants being slightly 
ho tter than implied by our r esults. The calibration 
bv lRarmrez fc M clcndea (|2005D is the only one that di- 
rectly provides color-temperat ure relations calibr a ted fo r 
giants. As already noted by ICasagrande et all ([2 010), 
the temperatures by iRamrrez fc Melendezl (|2005T ) seem 
to be systematically cooler than expected, and this is 
co nfirmed by ou r result s. Finally, the calibration given 
by IBruntt et al.l ((2012) overestimates temperatures at 
the cool end, which is again not surprising since their 
calibration was based on main-sequence stars only, and 
did not include corrections for lower surface gravities 
and different metallicities. Overall, we conclude that 
photometric estimates reproduce the measured temper- 
atures from interferometry well within the uncertain- 
ties, except for the giants where reddening is significant. 
We note that HD 173701 is the only star with sufficient 
Sloan photometry t o be included in the calibration by 
IPinsonneault et aTl ((20121 ). The SPS S temperature, cor- 
rected for metallicity as described in IPinsonneault et al.l 
(|2012f ). is 5364 ± 100 K, in good agreement to the de- 
termined values here. Finally, we note that the effective 
temperatures presented in this section do not influence 
the comparisons of the asteroseismic masses and radii 
calculated in the previous section (which were calculated 
using spectroscopic T e ff), since the dependence of Equa- 
tion ((4]) on T e ff is only small. 

4.3. Stellar Models 

Detailed modelling will be deferred to a future pa- 
per, but we present some first basic comparisons for the 
most interesting cases here. We use the publicly avail- 
able B aSTI stellar evolutionary tracks (|Pietrinferni et ahl 
I2004D with solar-scaled distribution of heavy elements 
(|Grevesse fc Noeli 119931 ) a nd a standard mass loss pa - 
rameter of r\ = 0.4 (see, e.g. iFusi-Pecci fc Ren zini 1976). 
The models do not include effects of diffusion or gravita- 
tional settling, and are calibrated to match the observed 
properties of the Sun with a mixing length parameter 
c*mlt = 1-913 and an initial chemical composition of 
(Y,Z) = (0.2734,0.0198). No convective-core overshoot- 
ing was included in the models presented here. Note that 
in the following we compare models to radii, masses and 
temperatures derived using our direct measurement of 
the angular diameter (see columns 7, 8 and 9 in Table [7]) 
and the spectroscopic metallicities. 

4.3.1. HD 182736 

The star with the best constrained fundamental prop- 
erties in our sample is the subgiant HD 182736, with rel- 
ative uncertainties in temperature, radius and mass of 
0.7%, 2.6% and 7.7%, respectively. The fact that the 
best observational result is achieved for the only sub- 
giant in our sample is not surprising: while the more 
distant red giants generally have well constrained diam- 
eters due to their larger size, they suffer from a large 
uncertainty in the parallaxes and effective temperatures 
due to their larger distance and significant reddening. 
On the other hand, main-sequence stars are generally 
too small to achieve a good precision on their measured 
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Fig. 10. — Radius versus effective temperature with the position 
of the subgiant HD 182736 shown as a red diamond. The solid 
line shows the BaSTI evolutionary model matching the metallic- 
ity from high-resolution spectroscopy and the mass determined in 
this study. Dashed-dotted and dashed-triple-dotted lines show the 
effect of varying the metallicity by 1 <r, while dotted and dashed 
lines show the same effect for varying the mass by 1 a. The deter- 
mined mass and metallicity for HD 182736 are M = 1.3 ± 0.1JW© 
and [Fe/H] = -0.15 ± 0.06. 



diameters. Subgiants land in the "sweet spot" between 
these regimes, with angular sizes big enough for a precise 
measurement with PAVO and distances close enough to 
have a well constrained Hipparcos parallax and negligible 
reddening. 

Figure [TU] shows a diagram of stellar radius versus ef- 
fective temperature with the position of HD 182736 ac- 
cording to the properties listed in Table [7] marked as a 
red diamond. The black solid line shows the evolution- 
ary track matching the determined mass and metallic- 
ity, calculated by quadratically interpolating the original 
BaSTI tracks. Dashed-dotted and dashed-triple-dotted 
lines show the effect of varying the metallicity by la, 
while dotted and dashed lines show the same effect for 
varying the mass by 1 a. The agreement between the 
models and our observations is excellent, with a match 
within 1 a for both radius and temperature. We em- 
phasize that no fitting is involved in this comparison - 
the mass, radius, temperature and metallicity are deter- 
mined independently from the evolutionary tracks. A 
more in-depth asteroseismic study using individual fre- 
quencies, in particular with resp ect to probing the cor e 
rotation rate using mixed modes (jDeheuvels et alJ (2012). 
combined with the results presented in this paper should 
yield powerful constraints for studying the structure and 
evolution of this evolved subgiant. 

4.3.2. HD 173701 

FigurelTTIshows the radius- T g diagram for HD 173701, 
a metal-rich main-sequence star with relatively well con- 
strained properties. In this case, the agreement between 
BaSTI models and observations is poor. The difference 
can be reconciled with a 3 — a difference in mass and 
metallicity, i.e. the star is more metal-rich and less mas- 
sive than implied from our observations. Indeed, the 
asteroseismic (but n ot m odel-independent) analyse s by 
iMathur et all ((20121) and iSilva Aguirre et al.l ((20121 ) im- 
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TABLE 7 

Fundamental properties of all Kepler and CoRoT stars in this study. 
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Vertical lines divide estimates based on spectroscopy (columns 3 
from estimates using the measured angular diameter (columns 7, 
reported for HD 176726 since our results suggest a measurement 
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Fig. 11. — Radius versus effective temperature with the position 
of the metal-rich main-sequence star HD 173701 shown as a red 
diamond. Lines compare BaSTI and CESAM evolutionary tracks 
with different masses and initial Helium fractions (see text). Note 
that each track starts at the zero-age main sequence (ZAMS). The 
determined mass and metallicity for HD 173701 are M = 1.05 ± 
O.O7M and [Fe/H] = +0.34 ± 0.06. 



ply a mass of 1.00 ± 0.01 M and 0.96 ± 0.04M© for 
HD 173701, respectively, which would significantly im- 
prove the agreement. We also note that the ~100K dif- 
ference to the spectroscopic T c g implies that the adopted 
metallicity may not be consistent with the interfero- 
metric T c g. However, as shown in Figure 111! even at 
the spectroscopic temperature of 5390 K the position of 
HD 173701 would still be slightly too cool for the mass 
determined from the interferometric radius and astero- 
seismic density. Additionally, adopting a lower T e s in 
the spectroscopic analysis would result in a lower metal- 
licity, and therefore enhance the disagreement between 
models and observations. 

A more interesting possibility is that the physical as- 
sumptions in the evolutionary models need to be adjusted 
to reproduce the properties of this star. To test this, we 
have computed additi onal tracks using the ID s tellar evo- 
lution code C ESAM (iMorel fe Lebretonl 120081). We use 
opacities from lFerguson et all (|2005| ) for the metal repar- 



and 4), asteroseismic scaling relations only (columns 5 and 6) 
8 and 9). No estimates based on asteroseismic constraints are 
error for this star (see text). 

tition bv lAsplund et all (120091). and NACRE nucle ar re- 
action rates are adapted from lAngulo et all ([1999). The 
models include diffusion and gravitational settling, and 
co nvection is described using the mixing length theory 
by iBohm-Vitensd ()1958[ ) with a solar calibrated value of 
1.88. We have computed two models with the spectro- 
scopically determined metallicity of [Fe/H] = 0.34 and a 
mass of 1.06-M©, once with solar-calibrated initial helium 
mass fraction Y — 0.2646, and once with Y = 0.2485, 
corresponding to the lower limit set by cosmological con- 
straints. Figure [TT] shows that changes in the initial 
chemical composition brings better agreement to our ob- 
servations. Similar changes can be invoke d by reduc- 
ing the mix in g- length parameter (see, e.g., IBasu et al.l 
l2010l 120121) . IWright et all (|200l list HD 173701 with 
a rotation period of 38 days and Ca H & K activ- 
ity of log(R' HK ) = —4.87. Both the slower rotation 
period and solar-like activity do not seem to be com- 
patible with a decreased convection efficiency (smaller 
mixing length parameter), which would be needed to 
bring the models in better agreement to our observa- 
tions. Additionally, a sub-solar helium mass fraction 
for HD173701 does not seem to be compatible with the 
roughl y linear helium-to-metal enrichment for metal-rich 
stars (Casagr ande et al]|2007|) . although the scatter in 
this relation is large and studies of the Hyades have con- 
firmed that stars can be depleted in hel ium and at the 
same time have a super-solar m etallicity (Lcb reton et al.l 
[200lHPinsonneault et al1[200l . 

An alternative explanation could be related to inad- 
equate modeling of stellar atmospheres for metal-rich 
stars. Systematics in these models would affect the 
bolometric flux and hence the determined effective tem- 
perature. Furthermore, systematic errors in the limb- 
darkening models for metal-rich stars would change the 
derived angular diameter, which influences the deter- 
mined radius, m ass and effec t ive te mperature. Detailed 
3-D models by IBigot et al.l (|2006[ ) for the metal-rich 
K-dwarf a Cen B showed less significant limb-darkening 
and hence slightly sm aller diameters compa red to sim- 
ple 1-D models, while [Chiavassa et al.l (|2010l ) found dif- 
ferences up to 3% for models of metal-poor giants. 
Such differences are expected to be enhanced in visibile 
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Fig. 12. — Same as Figure 1101 but for the red giant star 
HD 175955. The determined mass and metallicity for HD 175955 
are M = 1.1 ± 0.2M Q and [Fe/H] = +0.06 ± 0.15. 

wavelengths (such as the observat ions presented here) 
compared to infrared observat ions (|Allende Prieto et al.l 
l2002t lAufdenberg et all 120051 ). Furthermore, compar- 
isons of 1-D to 3-D models have also yielded higher fluxes 
for 3-D models, particularly at short wavelengths, which 
could lead to small inceases in the derived effective tem- 
perat ure (see, e.g.. lAufdenberg et al.l 1200 5; Cas agrandd 
2009). A higher effective temperature would bring bet- 
ter agreement to the evolutionary tracks and spectro- 
scopic estimates. More detailed modeling will be needed 
to confirm if refined estimates of limb- darkening, taking 
into account the non-solar metallicity for HD 173701, can 
explain the observed differences. 

4.3.3. HD 175955 

Figure [T2] presents a model comparison for HD 175955, 
a red giant with a well constrained angular diameter and 
the most precise Hipparcos parallax. Gravity mode pe- 
riod spacings measured using asteroseismology have been 
used to classify this sta r as a H-shell burnin g, ascending 
red giant branch star ([Bedding et al.1 120111 ). Figure LT21 
shows that the measured temperature of HD 175955 is 
slightly hotter than the position of the ascending RGB 
tracks, but overall in good agreement with its determined 
mass and metallicity. Similar to HD 182736, a combina- 
tion of the constraints presented here with detailed as- 
teroseismic studies (such as the measurement of mixed 
mode rot ational splittings to constrain the core rotation 
rate, see iBeck et al.l [20121 ) should allow a detailed the- 
oretical study of the internal structure and evolution of 
this star. 

4.3.4. Additional Notes 

We note that for a few stars the derived stellar proper- 
ties appear unphysical and are likely related to potential 
observational errors. For HD 175726, for example, the 
measured linear radius combined with the asteroseismic 
density implies a mass of 0.50 ± 0.03M Q , which seems 
incompatible with its measured radius, temperature and 
solar-metallicity. Using the spectroscopically determined 
metallicity and the radius and temperature from inter- 
ferometry, a comparison with BaSTI models indicates a 



mass of 1.07-M©, which would imply asteroseismic val- 
ues of Av ~ 142 fiHz and f max ~ 3300 /kHz. These val- 
ues are signif i cantly different than the results found by 
iMosser et al.l ()2009f ). The difference could be explained 
by an undetected companion causing a significant error 
in the parallax, or due to measurement errors in either 
the interferometric or the asteroseismic analysis. Unfor- 
tunately no CoRoT follow-up observations are planned 
for HD 175726, and hence a resolution of this discrepancy 
will have to await independent future observations. 

The red giant HD 181827, on the other hand, has a 
large mass which is difficult to reconcile with evolution- 
ary theory. Both the asteroseismic and interferometric 
constraints are solid, hence pointing to a potential prob- 
lem with the Hipparcos parallax. Indeed, HD 181827 has 
the largest fractional parallax uncertainty in our sample 
(10%), leading to a large uncertainty on the radius and 
hence mass. We note that HD 181827 has been astero- 
seismically identified as a secondary clump star (|Girardil 
1999; Beddin g et all 1201 ID . corresponding to a massive 
(> 2Mq) He-core burning red giant. Our result of a sig- 
nificantly higher mass for HD 181827 compared to typical 
red clump giants is hence qualitatively in agreement with 
its asteroseismically determined evolutionary state. 

5. CONCLUSIONS 

We have presented interferometrically measured angu- 
lar diameters of 10 stars for which asteroseismic con- 
straints are available from either the Kepler or CoRoT 
space telescopes. Combining these constraints with par- 
allaxes, spectroscopy and bolometric fluxes, we present 
a full set of near model-independent fundamental prop- 
erties for stars spanning in evolution from the main- 
sequence to the red clump. Our main conclusions from 
the derived properties are as follows: 

1. Our measured angular diameters show good agree- 
ment with the su rface-brightness relation by 
iKervella et"all (|2004l ) and the IRFM coupled with 
astero seismic constraints by iSilva Aguirre et all 
(|2012l ). with an overall residual scatter of 5%. 
Ou r results seem t o con firm that the relation 
by IKervella et al"1 (120041) and the method by 
ISilva Aguirre et all (|2012f ) are also reasonably ac- 
curate for red giants. 

2. A comparison of interferometric to asteroseismic 
radii calculated from scaling relations shows excel- 
lent agreement within the uncertainties. While the 
uncertainties for giants are large due to the uncer- 
tainties in the parallaxes, our results empirically 
prove that asteroseismic radii for unevolved stars 
using simple scaling relations are accurate to at 
least 4%. A test of the f max scaling relation also 
shows no systematic deviations as a function of evo- 
lutionary state within the observational uncertain- 
ties. 

3. A comparison of measured effective tempera- 
tures with estimates from modeling high-resolution 
spectra (mostly u sing the VWA method, see 
iBruntt et al.ll2010l) and f rom the photometric in- 
frared flux method (see ICasagran de et all I2010D 
shows good agreement with mean deviations of 
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-22 ± 32 K (with a scatter of 97 K) and -58 ± 31 K 
(with a scatter of 93 K), respectively, for stars be- 
tween T e s — 4600 — 6200 K. Some photometric cal- 
ibrations show slight systematic deviations for red 
giants, presumably due to the more significant in- 
fluence of reddening for these more distant stars. 

4. A first comparison of our results with evolutionary 
models shows very good agreement for the sub- 
giant HD 182736, while there appear to be some 
discrepancies for the metal-rich main-sequence star 
HD 173701. We speculate that these differences 
may be due to inadequate modeling of stellar at- 
mospheres or limb-darkening for metal-rich stars, 
but note that more detailed theoretical studies will 
be needed to confirm this result. 

While our study has demonstrated the potential of 
combining different constraints to test stellar model 
physics, it is clear that the overlap between the different 
techniques is still limited. This situation can be expected 
to be significantly improved with future projects such as 
the ground-base d network SONG (Ste llar Observations 
Network Group, IGrundahl et al.l I2006D . which will de- 
liver precise multi-site radial-velocity timeseries for as- 
teroseismology and exoplanet studies of nearby stars. 
On the oth er hand, the pl anned European space mis- 
sion Gaia (jPerrvmanl [2003) will provide accurate par- 
allaxes for stars down to V < 15, while potential up- 
grades of interferometers such as the CHARA Array with 
adaptive optics will push the sensitivity limits of inter- 
ferometric follow-up to fainter stars, therefore improv- 
ing the overlap with Kepler and future space-based mis- 
sion s such as T ESS ( Transiting Exoplanet Survey Satel- 
^ie. lRicker et al.ll2009D . The possibility of independently 
constraining radii, effective temperatures, masses and 
metallicities using asteroseismology, astrometry, interfer- 
ometry and spectroscopy for a large ensemble of stars to 
study stellar physics as well as to characterize potentially 
habitable exoplanets is clearly the next step for con- 
tinuing the exciting revolution induced by CoRoT and 
Kepler over the coming decades. 
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